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The Hikami-Larkin-Nagaoka (HLN) formula [Prog. Theor. Pliys. 63, 707 (1980)] describes the 
quantum corrections to the magnetoconductivity of a quasi-2D electron gas (quasi-2DEG) with 
parabolic dispersion. It predicts a crossover from weak localization to antilocalization as a function 
of the strength of scattering off spin-orbit impurities. Here, we derive the conductivity correction for 
massless Dirac fermions in 3D topological insulators (3DTIs) in the presence of spin-orbit impurities. 

We show that this correction is always positive and therefore we predict weak antilocalization for 
every value of the spin-orbit disorder. Furthermore, the correction to the diffusion constant is 


surprisingly linear in the strength of the impurity 
of experimental fits for the magnetoconductivity 
HLN formula. 

PACS numbers: 72.10.-d, 73.20.Fz, 73.43.Qt, 73.25.+i 

Introduction. The problem of the diffusion of the sur¬ 
face states of 3DTIs is a complex one due to a variety of 
competing phenomena. The most striking is the fact that 
these surface states are described by a Dirac Hamiltonian 
[ , 2], which gives rise to weak antilocalization (WAL) 
in the presence of scalar disorder [3-5]. The WAL cor¬ 
rection can be affected by the interaction of the surface 
states with the residual bulk states [6], the thickness of 
the film [7], or electron-electron interactions [8, 9] chang¬ 
ing its sign and turning it into weak localization (WL). At 
the same time, since 3DTIs have strong spin-orbit cou¬ 
pling, one expects spin-orbit coupled impurities to have 
a strong effect on transport, different however than in 
graphene where two valleys are present [10, 11]. Surpris¬ 
ingly this problem has received virtually no attention [12] 
and so far the Dirac nature of the states, that manifest 
itself in the angular dependence of the Green functions, 
has not been taken into account in studies of spin-orbit 
impurities [7]. The problem is even more complicated in 
a transverse magnetic field. 

The formula commonly used to fit the magnetocon¬ 
ductance experiments on 3DTIs [13—20] was derived by 
Hikami, Larkin and Nagaoka (HLN) [21]. The HLN 
formula, however, lacks important features relevant to 
3DTIs: it is derived for quasi-2DEGs with a parabolic 
electron dispersion (impurities are treated as three- 
dimensional objects), so it accounts neither for the Dirac 
nature of the surface states nor for their strictly two- 
dimensional character. Fig. 1 is the main result of this pa¬ 
per, and shows the different effects of spin-orbit scatter- 
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spin-orbit. Our results call for a reinterpretation 
of 3D TIs which have so far used the standard 



Figure 1: Comparison of the quantum correction to the con¬ 
ductivity as a function of the magnetic field for massless Dirac 
fermions (solid lines) and the Hikami-Larkin-Nagaoka formula 
for a quasi-2D system (dashed lines) for different values of 
the concentration of the spin-orbit impurities. Purple plots 
are for purely scalar disorder (A = 0), blue (green) plots for 
A = 0.2 (A = 1.25). The magnetic field is renormalized by 
Bo = e/(2nhvpTo), where to is the elastic scattering time 
in the absence of spin-orbit impurities and vf is the Fermi 
velocity. 


mg for the HLN formula, namely WL to WAL crossover 
with the strength of spin-orbit impurity scattering, and 
for Dirac fermions, where WAL appears regardless of 
the strength of spin-orbit impurities. We only observe 
a convergence of the two formulas for large values of the 
strength of the impurity spin-orbit scattering (when both 
equations are no longer valid) [28]. Moreover, the HLN 
formula gives the wrong value for the diffusion constant 
when the impurity spin-orbit coupling vanishes: it does 
not capture the fact that the absence of backscattering 
for Dirac fermions doubles the diffusion constant as com¬ 
pared to conventional electrons. 
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In this paper, we present the full calculation of coher¬ 
ent diffusive transport of Dirac fermions in the presence 
of both scalar and spin-orbit coupled impurities. We first 
show that, because of the winding of the spin around 
the Fermi surface, the Dirac nature of the surface states 
breaks the mirror symmetry around the :ry-plane (the 
disorder potential no longer commutes with the kinetic 
Hamiltonian), and thus allows for a correction to physi¬ 
cal quantities, such as the classical conductivity and the 
diffusion constant, which is linear in the strength of the 
disorder spin-orbit coupling as opposed to the quadratic 
dependence observed for a parabolic dispersion. More¬ 
over, we show that WAL appears for any strength of the 
spin-orbit coupled disorder for massless Dirac fermions, 
as opposed to the case of electrons with parabolic disper¬ 
sion, where a crossover from WL to WAL (no correction) 
is observed for quasi-2DEGs (strictly 2DEGs). 

Model. The Dirac Hamiltonian for the surface states 
of 3DTIs, including a random uncorrelated Gaussian dis¬ 
order potential comprising both scalar and spin-orbit im¬ 
purities, with impurity concentration rij, reads: 

n = hv F (k x a y - ky(T x ) + V(k, k'), (1) 


with V(k, k') = Uptp 1 + i\{k x k') ■ a the disorder 
potential. Here 1 is the 2x2 identity matrix in spin 
space, prr, = Ei e~ l ^ k ~ k k Rl with ‘I’ labelling the ran¬ 
dom locations of the impurities, and we have assumed a 
short-range impurity potential so that U is not a function 
of wave vector. The disorder correlator b(0 1 O') is 


(v(k, k')V(k\ k)^ = 7 o^Io+*AIi sin 7 +A 2 I 2 sin 2 7 ^, (2) 


We introduce the two characteristic times r and r* as 
y = 7 rp(£ , F ) 7 o(l + 4y) and £ = 7 rp(£ , F ) 7 0 A where 
p(E F ) = 2 o * s the density of states at the Fermi 
energy. We obtain for the imaginary part of the self en¬ 
ergy: —ImE(fc) = —— 1 LT ^ (cos 0& y — sin 0c x ). Near the 
Fermi surface the retarded and advanced Green functions 
G r / a = [E — (Ho ± E )] _1 take the form 


G R/A (fc) 


1 / 1 + cos 6a v — sin 0a x A 

2 \ E-hv F k±g- J’ 


( 5 ) 


where the corresponding mean free-time between two 
scattering events r e , derived through the Fermi golden 
rule, obeys the Matthiessen rule A = 1 | i = 

7 tp(-Ef) 7 o _|_ \ -|_ Due to the spin structure of both 

the Green function and the scattering potential, this elas¬ 
tic mean free-time shows an unusual linear dependence 
in the spin-orbit scattering strength A, as opposed to the 
case of non-relativistic electrons [21]. The mean-free time 
is no longer an even function of A, since the winding of the 
spin around the Fermi surface for Dirac fermions defines 
unequivocally the direction of the z-axis. 

Diffuson and renormalized current operator. The lad¬ 
der diagrams are responsible for the difference between 

the elastic scattering time r e and the transport time 
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r tr appearing in the diffusion constant D = L ' F J tr in a 
d-dimensional material. For Dirac fermions (graphene, 
3DTI surface states), in the presence of point-like scalar 
disorder scattering, a doubling of the transport time 
r tr = 2T e is observed. The diffuson structure factor r D 
obeys the Bethe-Salpeter equation (see Fig. 2) : 


where 7 = O' — 6 and 0 (O') is the polar angle of the wave 
vector k (k'). We have introduced 70 = n-JJ 2 to quantify 
the disorder strength, and the dimensionless parameter 
A = Afcp to describe the relative strength of the scalar and 
spin-orbit coupled disorder. We work in the limit A<C 1, 
meaning that the spin-orbit scattering length is much 
larger than the mean-free path. The I,; are: Io = l®!, 
Ii = er z 1 — 1 (g> cr z , and I 2 = a z ® <J Z . The Hamil¬ 
tonian (1) preserves time-reversal symmetry (TRS). Due 
to the non-commutativity of the Pauli matrices in the 
band and impurity Hamiltonians, terms linear in A af¬ 
fect charge and spin dynamics, as opposed to the case of 
spinless electrons where these linear terms are absent. 

The bare Green function reads: 


G° s Ak) = ~ 


1 ( l ss ' + cos 0a y , — sin Qa x 


E — hv F k ± *0+ 


( 3 ) 


The lifetime of a particle with wave vector k in the weak 
disorder limit is the imaginary part of the self energy E: 


r ^(M',g) = MM')+ (6) 

7<5 7 <5 

/ T^r^M'^G* {«')GA (A"-g)M0">0') > 

J (27 t) z 7 a f kS 

in which 0 (O') denotes the direction of the incoming (out¬ 
going) wave vector and Greek symbols are the spin in¬ 
dices. We will drop from now on the spin indices to 
simplify the equations. 

The Fourier decomposition of this diffuson structure 
factor as T D (0,0\q) = E n,m T n,m(q)e i{me '~ n9) allows a 
solution of the equation perturbatively in A. Up to second 
order, this diffuson structure factor acquires a non-trivial 
structure with angular components (e.g I\'j J , 7 ^ 0 ) that 
does not appear with isotropic scattering, even for Dirac 
fermions. This is due to the presence of higher harmonics 
in both the disorder correlator &(0,0'), and the Green 
functions (see supplementary material for details). 

The diffuson structure factor renormalizes the current 
operators j(k) (cf. Fig. 2) as: 


_ r xzy 

m = j ^2 b (o,o')G°(k'). 


(4) J(k')=j(k')+[G A (k)j(k)G R (k)T D (0,0\0). (7) 

Jk 
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Figure 2: Diagrammatic representation of the calculations. 
The first row represents the Bethe-Salpeter equation obeyed 
by the diffuson structure factor where the Greek symbols de¬ 
scribe spin indices. The second row shows the current renor¬ 
malization. The third row depicts the three Hikami boxes 
needed for the WAL calculation. 


We find explicitly for the renormalized current operator 
along the ^-direction 


qy2 

Jx(O') = -euF [2(1 - A + —)<J y -—- a + + h.c .], 

o o 

(8) 

where h.c. denotes the Hermitian conjugate and = 
a x ± ia v . For pure scalar disorder A 0 we recover the 
doubling of the current operator J x = 2 j x . Once again, 
the spin structure of the band and disorder Hamiltonians 
is responsible for a non-trivial angular dependence of this 
renormalized current operator, no longer proportional to 
the original current operator j x = —ev F a v . 

Longitudinal conductivity and diffusion constant. It is 
now possible to compute the longitudinal conductivity 
through the Kubo formula as : 


0.Dr _ 


2ttLI 


ReTr 


J x (k)G R {k)j x {k)G A {k ) , (9) 


where Tr denotes a sum over both spins and wave vector 
k. and fl is the volume. Up to second order in A, we find : 


= e 2 p(E F )vlT e ^1 - A + . (10) 


The diffusion constant D = v F r e ^1 — A + shows a 
linear dependence on the spin-orbit scattering strength, 
as distinct from the case of non relativistic electrons. As 
this diffusion constant is a crucial parameter in weak an¬ 
tilocalization, we expect the behavior of Dirac fermions 
to be different from the usual HLN formula [21]. 

Weak antilocalization correction. The weak antilocal¬ 
ization correction requires the three different contribu¬ 
tions pictured in Fig. 2, as is usual for Dirac fermions [3]. 


The Cooperon structure factor r c (maximally crossed di¬ 
agrams) accounts for the quantum interferences of closed 
paths during diffusion. Because of TRS, it can be ob¬ 
tained from the diffuson structure factor by twisting the 
retarded or advanced branch. It obeys: 

T%{6ff',q) = b c af} {e,0')+ ( 11 ) 

7<5 7<5 

[ c a u(eff"^G R ffV')G A ffq-V')b%{6"ff') , 

J \ Z7r )~ 7* kS 

where b c = 7 0 ^I 0 + iAIisiny — A 2 I 2 sin 2 7 ^ and R = 

a z i0 1 + 1 g> o z . We also expand r c (0, O') = 
keeping only the 9 terms in 1/g 2 

up to second order in A (the terms with n, m £ 0,±1) 
(see supplementary material). 

It is then possible to calculate the three Hikami boxes 
pictured in Fig. 2 for each of these nine terms (full ex¬ 
pressions in supplement). The contribution of each r([ rn 
mode can be collected in three different groups, depend¬ 
ing on the value i = \n\ + |m| = 0,1, 2. The respective 

• •• 1 2 \ ^ 

weights Wi of these contributions are wq = ^ ~ A + £ 5 -, 

Wi = j —and W 2 = -g-. Summing these 9 contribu¬ 
tions, we obtain the WAL correction : 


5f7DF = ^ ln (^) • (12) 

This expression should be compared with the formula 
for non-relativistic electrons for strictly 2D and quasi-2D 
systems : 


c (2D) _ 

®°”hln — 


— In 


1 + A 2 /2 


5 a 


(q-2D) 

HLN 


7 ih 


£ + A V 2 

e 2 

In 


7T h 

W + RT 




1 + ^L T A 


(13) 

These three formulas are plotted in Fig. 3 as a func¬ 
tion of A. We have renormalized these corrections by 
5a 0 = hi ( j, where r e depends on the model and 


is a function of A. We have set the ratio — where rn 

jo . u 

is the value of r e in the absence of spin-orbit scattering 
To = 7rp (J F ) 7o to be equal to 10 in agreement with what 
is measured experimentally[22-24]. We observe that the 
Dirac fermions remain in the same symmetry class (sym- 
plectic, with WAL), whereas the HLN formula shows a 
crossover from the orthogonal symmetry class (WL) to 
either no correction for strictly 2DEG, or WAL for quasi- 


2DEG. 


Due to the renormalization of the scattering time by 
spin-orbit impurities, one could also interpret Eq. (12), as 
increasing with A if one renormalizes conductivity correc¬ 
tions by the scattering time in the absence of spin-orbit 
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Figure 3: Comparison of the quantum correction to conduc¬ 
tivity as a function of the spin-orbit scattering strength A for 
massless Dirac fermions (solid line), Hikami-Larkin-Nagaoka 
formula for strictly 2D system (dashed line) and quasi-2D 
system (dotted line). The conductivities are given in units of 
Sag = In To emphasize the dependence of WAL on 

the spin-orbit impurity strength for massless Dirac fermions, 
the inset shows WAL in units independent of A i.e. renormal¬ 
ized by ^ In (to is the scattering time in the absence of 

spin-orbit impurities). The ratio /to varies over two decades 
from 5 (purple), 10 (blue), 50 (green), 100 (orange) to 500 
(red). 


impurities r 0 (namely normalized by lnr^/ro) : 


<5fD.F 


2nh 



+ A + 0(A 2 ) 


(14) 


The inset to Fig. 3 shows the linear dependence of the 
normalized correction to the conductivity as a function 
of A for massless Dirac fermions. As A = Aone can 
experimentally probe this linear dependence by varying 
the Fermi wave vector using an electrostatic gate, for a 
constant strength of spin-orbit impurity scattering. 

WAL correction as a function of the magnetic field. 
The well-known F1LN formula [21] describes the quan¬ 
tum correction to conductivity as a function of an ap¬ 
plied magnetic field for non-relativistic electrons in pres¬ 
ence of both scalar and spin-orbit impurities, where the 
only relevant parameter is the diffusion constant. In 
contrast, in our problem, each mode obeys a dif¬ 

fusion equation with a diffusion constant Di depending 
on i = |n| + |m|. Namely, these diffusion constants are 
D 0 = ti|r e /(l-|-A 2 / 2 ), D 1 = 2u|r e /A and D 2 =2u|r e /A 2 . 

Inserting the magnetic field through the Peierls sub¬ 
stitution maps this diffusion equation onto a Schrodinger 
equation with the equivalence Di h/2m and e <->• 2 e 
[25-27]. It follows that the introduction of the magnetic 
field reduces the contribution of each mode to the WAL 
correction by a factor : 


9i(B) = 




(* + ?)-*(* + ¥) 

In a 


(15) 


where we define the characteristic fields Bf = and 
Bf = Collecting all the modes we obtain : 


5a D . F .(B) = -^yln — 'Y\a i w i g i {B) , 

7 Til T 4 / 


T _± 

To 


(16) 


»=l 


where the mode weights at are ai = 1, and 02 = 03 = 4, 
with the w,s defined before Eq. (12). The magnetocon¬ 
ductivity corrections for massless Dirac fermions are pre¬ 
sented in Fig. 1. 

Discussion. Our results show many differences from 
the HLN formula widely used to fit magnetotransport 
experiments of 3DTI surface states. The first is that the 
winding of the spin around the Fermi surface breaks the 
mirror symmetry around the xy- plane, so it is now pos¬ 
sible to obtain a linear dependence of measurable quan¬ 
tities on A, the strength of the spin-orbit disorder. Such 
a linear dependence is observed in the mean-free time r e , 
the longitudinal Drude conductivity cr xx , and the diffu¬ 
sion constant D. 

A second type of difference emerges from the 
anisotropy of the Green functions for massless Dirac 
model. This anisotropy together with the anisotropy 
coming from the spin-orbit impurities leads to 9 different 
Fourier modes in the Cooperon to the second order in A as 
opposed to only one mode for the HLN model. Moreover, 
our expansion in the spin-orbit impurity strength shows 
more explicitly the fact that this calculation is perturba¬ 
tive in A and should be restricted to small values of the 
perturbative parameter as the odd powers of the series 
expansion contribute negatively to the conductivity. 

In general symmetry terms, the massless Dirac fermion 
model stays in the symplectic class for all values of the 
impurity spin-orbit coupling, as the square of the time 
reversal operator 0 2 = —1. This explains why WAL is 
always observed for Dirac fermions. In the HLN formula, 
the introduction of the impurity spin-orbit coupling is re¬ 
sponsible for a crossover from the orthogonal class 0 2 = 
1 when A = 0 to the symplectic (’’pseudo-unitary” [29]) 
class when A —» 00 for quasi-2DEGs (strictly 2DEGs). 

Conclusions. We have derived the magneto conductiv¬ 
ity corrections for the surface states of 3D TI in the pres¬ 
ence of scalar and spin-orbit impurity disorder. This is 
expected to be directly relevant to the experimental anal¬ 
ysis of these materials. We showed the profound differ¬ 
ence between the HLN model for parabolic dispersion 
and the massless Dirac fermion model. For the latter we 
predict WAL in the presence as well as in the absence 
of spin-orbit impurity scattering, while the diffusion con¬ 
stant and the longitudinal conductivity are renormalized 
to the first order in the spin-orbit scattering strength. 
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Supplementary material 


Solution of the Bethe-Salpeter equation for the diffuson 


The calculation of the classical conductivity requires the current operator renormalization, through the insertion of 
the diffuson, also known as the ladder diagram. To obtain this diffuson, we need to solve the Bethe-Salpeter equation : 

/ dk" _ _ 

j^^x(0,e'',q)G^(k")Gi x (k'' - q)Kp, KS (6",6') , (17) 

where 9 (O') denotes the direction of the incoming (outgoing) wave vector and Greek symbols describe the spin indices. 
Our first step is to expand all the quantities in Fourier series in order to remove the angular dependence. We write : 


r^ s (0,0',q) = J2 T ?n, m)a ^s(<f)e l{me '- nB) 

n,m 

W(M') = J2 b (n,m) a ^Se i{m8 '- ne) , 


(18) 

(19) 


as b only depends on O' — 0 , we will use b^ nrrF) = b^Sn^m where 6 is the Kronecker symbol. Looking at each mode in 
the Fourier expansion, we obtain a system of coupled equations : 

■^(n,ra)a:/3,7<5 (5) = b { n)a/3,7<5^n,ra 

k 


where Pf k) ^ XK (q) is defined as : 


pD 

1 (k)fiiy,\K, 


(q) 


J ^e ifee "G« (k")G^(k" - q) (21) 

_1_ f do" (1 +cos Q"gV - sin 0"a x )^ ® (1 - cos B" a* - sin 0" <t x )xk cik9 » ^ 

2 7o f 1 + A + A 2 /2p 2n 1 + iv F T e (q x cos 0” + q v sin 0”) 


It is important to highlight the transposition made in the advanced Green function in order to write this new 
Bethe-Salpeter equation as a matrix product, which allow to drop the spin indices from now on. As we are mainly 
interested in the small q, we perform a Taylor expansion of Eq. 22 around q = 0 : 

Pfk)(o) = Pfk) - ™FT e (q+Pf k „ 1 ) + q-Pf k+ 1 )) - vlrl (q 2 + Pf k _ 2) + 2q+q_pf k) + q 2 _Pf k+2) ) . (23) 

In the absence of spin-orbit impurities, only bo is non-zero, and the only contribution to the ladder diagram is through 
Pq. When the spin-orbit impurities are present and taken into account, only 5 b n are non-zero, and 5 P n contributes 
to second order in A, for n = 0, ±1, ±2. 

Now to solve the coupled equations system, we use the fact that the disorder correlator only shows a small number 
of harmonics ( b n = 0 for n ^ 0,±1,±2), and the fact that the Eq. 20 states that F^ m oc b m to ensure that not all 
the modes will contribute. Moreover as b n oc A n l, a series expansion in powers of A is a natural choice. We write all 
the quantities as X D = X^X a , and solve the system iteratively for every order of A : 

r { n a) m (q)=b^6 n , m + £ r ( :i + M P i a2 \q)bt 3) ■ (24) 

k,a 1 +a 2 +a 3 =a 


To zeroth order in A (so in the absence of the spin-orbit scattering), only the will contribute (it corresponds to 
the case of Dirac fermions in the presence of scalar disorder that has already been studied many times). Then we look 
at the first order in A, and we calculate T±i : ±i, r±i ; o and also the first order contribution to To^- We have calculated 
this diffuson modes up to the second order in A, and obtained 15 different modes T n jn with n,m = 0, ±1 ± 2 (25 
modes minus the 10 components F-j-i, r± 2 ,±i and F_|-( 2 ,- 2 ) which are of higher order in A). For example, we obtain 
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in the up-down basis along the z-direction : 

( 0 0 0 0 ^ 

4, 0 0 ii 

X' 0 0 X' 

\ 0 0 0 0 / 

a x 0 t + a x 0 a z — io y 0 3L — 0a 2 1 0 a x + il 0 a y — a z 0 cr* — ia z 0 <r y 


-i,o — 


(25) 


= X' 

+^i 

with X' = ( 2 \ 2 - 

2 \VpT^q 2 16 


4 4 

o- x (g) 1 — g> cr z + ia y ® 1 — ia y (g> ct z | t ® o x - it ® <j v + a z ® cr x ~ ia z ® 
4 1 4 

and A = ~ i H~ 2U (1 — A). 


(26) 


Expression for the Cooperon modes and their contribution to WAL 

We can solve the Bethe-Salpeter equation of the Cooperon with the same technique, or we can use the symmetry 
between the diffuson and the Cooperon (the Cooperon correspond to a diffuson where the advanced Green function is 
time-reversed). However, we know that the main contribution to the conductivity will be given by the terms in 1/q 2 , 
so we keep only the modes with such terms, the other ones being negligible for diffusion at long distances. We obtain 
the 9 modes : 


r 


c 

(0,0)a/3,7<5 


r 


c 

(±l,=Fl)a/3,7<5 


r 


c 

(±l,0)a/3,7<5 


1 (0,±l)a/3,7<5 


F c 

1 (=t=l,=bl)a/3,7<5 


1 + A 2 /4 lap ® 17-5 - <0 ® <%s - <8 ® - <P ® g 7 s . 

VpT^q 2 4 

A 2 70 <0 ® a *6 ± ia l0 ® ± l <0 ® <**S ~ *1(1 ® tfs . 

2 VpT 2 q 2 4 

-iAyo ( To^p ® 1 7 « - a x p ® a z 5 - ia v a0 ® 1 74 =F ^ ® ** 4 


2vlr 2 q 2 


ila/3 ® <T* S ± la/3 ® C7* 4 + (7^ ® 0^ 4 ± *<7^ ® <7* 4 ' 


— *A 7 o / ±^S/3 ® + <0 ® 0-75 - i(jV a0 ® *7* ? ia lp ® ^75 


2^|r e V 


«l a /3 ® ^ T la/3 ® ^ ~ ^ ® ^ ± ^ ® ^' 


A 2 7 o la/3 ® 1 7 <5 ± la/3 ® ^5 ± <0 <g> 1 7<5 + <7^ <g> C7 * 4 
2 VpT 2 q 2 4 


(27) 

(28) 


(29) 


(30) 

(31) 


Each one of these 9 modes contributes to the weak anti-localization when included in the three Hikami boxes (the 
bare one and the two dressed ones) pictured in Fig. 2 of the main text. For each mode, these three contributions can 
be written in the form 5a = w f ' h In The following table gives the value of the weight w for each of the mode, 

and for each of the Hikami boxes. 


Mode 

Bare H.B. 

Dressed H.B. (x 2) 

Total iU| n |+| m | 

•I 

o 

o 

1 — 3A + (25/4)A 2 

— (1/2) + 2A — (19/4)A 2 

(1/2) — A + (3/2)A 2 

r n,m (|rc| + \m\ = 1) 

(1/2)A — (3/2)A 2 

— (1/4)A + A 2 

(1/4)A — (1/2)A 2 

r„, m (n + m 1=2) 

(1/4)A 2 

(1/8)A 2 

(1/8)A 2 

Total 

1 - A + (5/4)A 2 

(1/2) + A — (5/4)A 2 

1/2 


This result is plotted in Fig. 3 of the main text, and shows that it is necessary to take into accounts all the Fourier 


modes of the Cooperon to obtain the conductivity correction do- 



characteristic of the symplectic class. 





Figure 4: Crossover from WL to WAL obtained by the Hikami-Larkin-Nagaoka formula for a quasi-2DEG for different values 
of the ratio t^/to over two decades from 5 (purple), 10 (blue), 50 (green), 100 (orange) to 500 (red). The conductivities are 

given in units of 8ao = ^ In (j^)- 


Range of validity of the HLN formula 

Fig. 4 shows the crossover from the weak localization to the weak antilocalization for quasi-2DEG with parabolic 
dispersion when the spin-orbit impurity concentration increases. For metals, where the ratio r^/re is very large (of 
the order of 1000 [27]) this crossover occurs for a value of A small enough that a perturbative treatment is possible. 
However, for the parameters experimentally relevant for 3DTIs, with a smaller ratio r^/r e around 10 [22-24], this 
crossover occurs for values of A of the order of the unity, which is beyond the range of validity of the HLN derivation. 
We have plotted in Fig. 4 the crossover from WL to WAL for different values of the ratio r e to highlight that the 
formula derived by Hikami et al. is not enough to explain the WAL correction observed experimentally in 3DTIs. 





